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Abstract 

The existence of intimate relation between generalized statistics and supersymmetry 
is established by observation of hidden supersymmetric structure in pure parabosonic 
systems. This structure is characterized generally by a nonlinear superalgebra. The 
nonlinear supersymmetry of parabosonic systems may be realized, in turn, by modifying 
appropriately the usual supersymmetric quantum mechanics. The relation of nonlinear 
parabosonic supersymmetry to the Calogero-like models with exchange interaction and 
to the spin chain models with inverse-square interaction is pointed out. 
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1 Introduction 



One of the fundamental and seminal approaches in theoretical physics is a search for sym- 
metries. When new symmetry is revealed, this leads not only to better understanding of a 
system but sometimes gives rise to establishing unexpected and important relations between 
different theories and lines of research. 

Generalized statistics first was introduced in the form of parastatistics as an exotic pos- 
sibility extending the Bose and Fermi statistics [jTJ-PJ and for the long period of time the 
interest to it was rather academic. Nowadays it finds applications in the physics of the 
quantum Hall effect and (probably) it is relevant to high temperature superconductivity 0. 
Super symmetry, instead, unifies Bose and Fermi statistics 0-[(| and its development lead to 
the construction of field and string theories with exceptional properties JTOl 11], that trans- 



formed the same idea of supersymmetry in one of the cornerstones of modern theoretical 
physics. Though supersymmetry and generalized statistics may be unified in the form of 
parasupersymmetry JT2 



TBI 



nevertheless, by the construction, the two concepts seem to be 
independent. 

In this paper we establish the existence of intimate relation between generalized statistics 
and supersymmetry by observation of hidden supersymmetric structure in pure parabosonic 
systems. This structure is characterized generally by a nonlinear superalgebra. The nonlin- 
ear supersymmetry of parabosonic systems may be realized, in turn, by modifying the usual 
supersymmetric quantum mechanics. More specifically, we investigate hidden supersymme- 
tries of the single-mode parabosonic systems described by the Hamiltonians of the simplest 
quadratic in parabosonic creation-annihilation operators form and show that such simple 
systems possess a reach set of N = 1 supersymmetries. Depending on the order of a parabo- 
son, supersymmetry is realized either linearly in the unbroken or spontaneously broken phase 
II H, or in the form of nonlinear symmetry generalizing for the supersymmetric case 

the known nonlinear symmetries of the system of harmonic bosonic oscillators with rational 
ratio of frequencies f2(J. The basic tool with which we analyze the hidden supersymmetries 
of pure parabosonic systems is the so called .R-deformed Heisenberg algebra (RDHA) ||21||- 
25| , introduced inexplicitly by Wigner under generalization of the bosonic oscillator that 
lead subsequently to the concept of parafields and parastatistics [IIJ, 0, Hj,p6f-p8f. Due to 



the RDHA structure underlying parabosons, as we shall see, the systems under investiga- 
tion turn out to be related to inverse- square interacting particle systems ||, |3D|, which 
are important models of many-body systems because of their exact solvability and intimate 
connection to spin chain systems [[H]] and 2-dimensional Yang-Mills theories |3"2"| |. Recently, 
this special deformation of Heisenberg algebra has found applications for the description 
of 3D f33[-|3(| and 2D |37| anyons, 3D supersymmetry and higher-spin gauge interac- 
tions PSJ . It lies also in the basis of the construction of the minimally bosonized version of 



supersymmetric quantum mechanics 34, 39, 40 



One of the basic ingredients of supersymmetric quantum mechanics 0, || is the Z2- 
grading operator T, T 2 = 1, which classifies all the operators into even (bosonic, B) and odd 
(fermionic, F) subsets according to the relations [T, B] = {T, F} = 0. The main idea of the 
minimally bosonized supersymmetric quantum mechanics is to take the reflection (parity) 
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operator R, Rty{x) = x), as a grading operator. With this, the operators 



Q l = tRQ 2j Q 2 = -^(-^ + W{x)r\, (1.1) 
<h + W\x) - W'(x)r) , (1.2) 




containing an odd function W(x) = —W(—x), are the odd supercharges, {R, Qi} = 0, 
i = 1,2, and the even Hamiltonian, [R, H] = 0, satisfying the linear superalgebra of Witten's 
supersymmetric quantum mechanics: 

{Q + ,Q-} = H, Q ±2 = 0, [H,Q ± )=0, (1.3) 

where Q ± = \{Q 1 =F«Q 2 )- 

With the choice of the superpotential W(x) = — j-, v e R, the operator \f2Q 2 coincides 
with the Yang-Dunkl operator (41| 



D " = - i (s-s JI )' (L4) 

related to the Calogero model 0, p2]| , where i? plays the role of the exchange operator. With 
the extended differential operator D u , one can construct the analogs of bosonic creation- 
annihilation operators, 

a ± = —7={x =F iD u ). (1.5) 
v2 

These operators together with the reflection operator R form the -R-deformed Heisenberg 
algebra, 

[ a -,a + ] = 1 + pR, {R,a ± } = 0, R 2 = 1, (1.6) 

which possesses unitary infinite-dimensional representations for v > — 1 and at the integer 
values of the deformation parameter, v — p — 1, p — 1, 2, . . ., is directly related to parabosons 
of order p 0, 23, 2f|. On the other hand, at v — —(2p + 1) the algebra (|1.6|) has finite- 



dimensional representations corresponding to the deformed parafermions of order 2p (see 
Appendix). 

Earlier it was observed [TJ9] that the systems given by the Hamiltonians 



H e = ^{a + ,a-}-^eR[a-,a + ], e = +, -, (1.7) 



are supersymmetric: they are described by Eqs. ( |1.1|) , ( |1.2|) with W = ex — j-. These 



systems possess the following property: in the case e = + the spectrum of the Hamiltonian 
( P..7D does not depend on the value of the deformation parameter u > —1, i.e. Eq. (|1.7|) gives 



the isospectral family of supersymmetric systems being in the phase of exact supersymmetry 
{Eq = 0). On the other hand, at e = — the Hamiltonian (|1.7|) describes the set of systems in 
the phase of spontaneously broken supersymmetry with the scale of supersymmetry breaking 
Eq = 1 + v > governed by the value of the deformation parameter v > — 1 . 
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Here we first show that at v = 1 the unbroken and spontaneously broken supersymme- 
tries of the system (|1.7|) are the linear supersymmetries of the parabosonic systems described 
by the Hamiltonians of the simplest form H n = a + a~ and H a = a~a + . Then we shall find 
that at v — 2k + 1, k — 1, 2, the same quadratic Hamiltonians supply us with a super- 
symmetry characterized by nonlinear superalgebra. The peculiarity of linear and nonlinear 
supersymmetries of pure parabosonic systems, as we shall see, is encoded in the nature of 
their supercharges: being realized in terms of parabosonic creation-annihilation operators 
a ± , they are essentially nonlinear. We shall observe that at the values of deformation pa- 
rameter v = 2k corresponding to the parabosons of order p = 2k + 1, the spectra of the 
Hamiltonians H n and H a have no supersymmetry but possess the interesting structure: they 
reveal a finite number of 'holes' in comparison with bosonic spectra [y = 0). This 'hole num- 
ber' is correlated with the parastatistics' order. It will be shown that with the single-mode 
supersymmetric parabosonic systems one can associate the systems of two identical fermions, 
which are the simplest two-particle spin chain models with inverse- square interaction. We 
shall also discuss the realization of parabosonic supersymmetries by modifying the classical 
analog of usual supersymmetric quantum mechanics. 

The paper is organized as follows. Section 2 is devoted to investigation of the supersym- 
metries of the parabosonic systems given by the normal ordered Hamiltonian H n = a + a~ . 
Section 3 deals with the systems described by the Hamiltonian H a = a~a + . Here, in partic- 
ular, we find the exotic parabosonic system with the spectrum of usual superoscillator ||2| 



but characterized by the nonlinear superalgebra. Section 4 concerns the realization of the 
parabosonic supersymmetries in the system of two identical fermions. In Section 5 we show 
how the nonlinear supersymmetry of pure parabosonic systems may be reproduced by appro- 
priate modification of the classical analog of the usual supersymmetric quantum mechanics. 
In concluding Section 6 we discuss the revealed supersymmetries from the view-point of gen- 
eralized deformed oscillator ||43|| , Calogero-like models with exchange interaction 0, |30, ^4 



reduced parasupersymmetry, and nonlinear finite ly-symmetries and give some list of 
open problems for further investigation. Appendix contains the necessary information on 
irreducible representations of RDHA and on the relationship of RDHA to parabosons and 
parafermions. 



2 The systems with Hamiltonian H n = a + a 

Let us choose the superpotential in Eqs. ( |1 . 1|) , (|1.2|) in the form 

W = x-— . (2.1) 

In this case the Hamiltonian (|1.2j ) can be represented as the normal ordered product 

H n = a + a~ (2.2) 

of the operators ( |1.5|) corresponding to the particular value of the deformation parameter: 
v = 1. This form of the Hamiltonian can also be obtained from Eq. ( |1.7| ) by putting in 
it e = + and v — 1. The supersymmetric Hamiltonian (|2.2| ) may be written equivalently 
as H n = N + n_, where iV is a number operator, [N, a^] = ia*, and n_ = |(1 — R) is 
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a projector (see Appendix). From here the spectrum of the normal ordered Hamiltonian 
( p. 2D with ai^ corresponding to the case v = 1 of algebra (|1.6| ) can be immediately obtained: 
H n \l) = Ei\l), E = 0, E21-1 = E21 — 21, I = 1,2,.... So, the bosonized supersymmetric 
system characterized by superpotential (|2.1| ) is the single-mode parabosonic (of order p — 2) 
system with unbroken supersymmetry. 

Let us take the normal form for the Hamiltonian ( |2.2| ) in general case of unitary represen- 
tations of RDHA [y > —1) and investigate its spectrum. In terms of the number operator 
we have 

H n = N + uU- 

The even states \2l) are stable eigenstates of H n : their eigenvalues do not depend on the 
deformation parameter, 

whereas the eigenvalues corresponding to the odd eigenstates |2/ + 1) are subject to the linear 
shift, 

E» l+l =E^ + v = 2l + l + v. 

There are two special cases given by v = 2k and v = 2k + 1 with k = 1,2,.... For v = 2k 
the spectrum of the system (|2.2|) contains k 'holes' in comparison with the non-deformed 
(bosonic) case v = 0: the first k odd energy levels with E = 2h — 1, h = 1, . . . , k, are absent 
from it, whereas the rest of the spectrum is the same as for the bosonic system (see Fig. la). 

On the other hand, the odd values of the deformation parameter, v = 2k + 1, supply us 
with supersymmetry. As we have seen, at k = the Hamiltonian (|2.2| ) describes the system 
in the phase of unbroken supersymmetry with one singlet state corresponding to E — 
and all other states to be paired in supersymmetric doublets. For the parabosons of order 
p = 2(k + 1), k — 1,2,..., one gets the nonstandard supersymmetry. It is specified by the 
presence of k + 1 singlet states |2Z), / = 0, k, with energies E 2 i = 0, 2, ... , 2k, whereas 
other states are paired in supersymmetric doublets with £? 2 (/-/s)-i = E 2 i = 21, I — k + 1, . . . 
(see Fig. lb). 

One notes that the case of the nonstandard N = 1 supersymmetry with several sin- 
glet states may be treated as a reduced parasupersymmetry. Indeed, let us consider the 
parasupersymmetric system described by the Hamiltonian 

H = b + b- + J 3 +j. (2.3) 

Here b are the bosonic creation-annihilation operators and J3 is one of the su(2) generators 

taken in spin-j representation, J 2 = j(j + 1), with j = (k + l)/2, k = 1, 2, The reduction 

of this spin-j parasupersymmetric system by two additional linear constraints 

(J 3 -j)l*> = 0> (J 3 +m)=0 (2.4) 

(or by one quadratic constraint (Jf — j 2 )\^) = 0) results in the spectra corresponding to the 
nonlinear supersymmetry with k + 1 singlet states (see Fig. lc). 

Let us investigate in more detail the nature of the observed nonstandard supersymmetry. 
First, one constructs the odd generators 

Q + = a + n_, Q-=a-U+, (2.5) 
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Figure 1: a) Spectrum of the system ( |2.2| ) at v — 4; b) Spectrum of the system ( j2.2| ) at 
v = 5; c) Spectrum of the reduced parasupersymmetric system ( |2.3|) , fl2.4|) for j = 3/2 with 
the states j'3 = ±1/2 eliminated by constraints ( |2.4|) . 



with ri-t- = ^(1 ± i? ). Together with the Hamiltonian and grading operator R they satisfy 
the following relations: 



Q 



±2 



0, {R,Q ± } = 0, [R,H n ) = 0, 



{Q + ,Q-} = H n + k(R-l), [H n ,Q d 



T2kQ z 



±2Q" 



R 2 



1. 



[R,Q±\ 

This can be understood as a central extension of the supersymmetric algebra ( |1.3| ) with the 
grading operator R nontrivially included in the set of even generators. Only for k = the 
odd generators Q ± commute with H n and are the odd integrals of motion (supercharges). 
In this case they coincide with the operators Q ± = \{Q\ =F iQ-z), where Qi,2 are constructed 
according to ( |1 . 1| ) with W = x — g - . Therefore, we arrive at the natural question: what are 
the supersymmetric generators being odd integrals of motion which mutually transform the 
supersymmetric doublet states and what kind of superalgebra do they form? Such mutually 
conjugate odd operators satisfying the relations {R, = are 



q+ = (a + ) 2fc+1 n_, q- = (a-) 2fe+1 n H 



(2.6) 



Together with the Hamiltonian H n they form the following nonlinear superalgebra: 

Q ±2 = 0, {Q + ,Q~} = H k (H n ), 1^,0^=0, (2.7) 

where 



n k {H n ) = (H n - 2k)(H n -2k + 2)...H n ...(H n + 2k - 2){H n 
This algebra follows immediately from the relations 

Q + Q- = n k (H n )U + , Q-Q+ = H k (H n )U_. 



2k). 
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The singlet states |0)j = \2i), i = 0, 1, . . . , 2k, are annihilated by the odd generators, Q ± \0) i = 
0, but only the first state |0) = |0) is the zero energy ground state: H n \0) = 0. Thus, in the 
case v = p — 1 = 2k + 1 the supersymmetry of the single-mode parabosonic system ( |2.2| ) is 
characterized by the presence of k + 1 supersymmetric singlet states and for parabosons of 
order p — 4, 6, . . ., unlike the case of p = 2, the corresponding superalgebra is nonlinear. 

Realizing the creation-annihilation parabosonic operators via Eqs. (|1.4fj , (|1.5|) , we repre- 
sent the Hamiltonians H n = a + a~ and supercharges (|2.(j|) in the form 



Hn = ^(-^+^ + ^-1+^7^-1)^), (2- 




2fc+l 



with p = 2k + 1. The system given by the Hamiltonian (|2.8|) can be treated as a 2-particle 
Calogero-like model with exchange interaction, where x has a sense of a relative coordinate 
and i? has to be understood as an exchange operator. The systems of such type were 



considered by Polychronakos ||, [ft], |44|. Therefore, at odd values of the parameter v the 
class of Calogero-like systems ( |2.8| ) possesses a hidden supersymmetry which at v — 1 is the 
N — 1 linear supersymmetry in the unbroken phase, whereas at v — 2k + 1, k — 1, 2, 
the supersymmetry is characterized by the supercharges being differential operators of order 
2k + 1 satisfying nonlinear superalgebra 



3 The systems with Hamiltonian H a = a a + 

Let us choose now the superpotential in Eqs. (|1 . 1|) , (|1 . 2| ) in the form 

W = -x - — . 

2x 

In this case the Hamiltonian can be represented equivalently as the antinormal product 

H a = a~a + (3.1) 

of the creation-annihilation operators of RDHA corresponding to the case of parabosons 
of order 2 [y — 1). This antinormal form may be obtained also from the supersymmetric 
Hamiltonian ( |1.7| ) by putting in it e = — , v — 1. The corresponding mutually conjugate 
operators Q ± = \{Q\ =F iQ-z) take now the form 

Q + = a + n+, QT = a~n_ (3.2) 



instead of supercharges ( [2.5|) of the system (|2.2|) at v = 1. Representing the Hamiltonian ( |3.1|) 
in the equivalent form H a = A^+l+n + in terms of the number operator and projector n + , we 
conclude that it describes the system in the phase of spontaneously broken supersymmetry 
characterized by the spectrum E\ = 2 [Z/2] + 2, / = 0, 1, . . . , where [.] is the integer part. 

Let us follow the same line as in the previous section and treat the creation-annihilation 
operators in the Hamiltonian ( ^.1| ) as the generators of RDHA for the case of generalized 
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parabosons [y > — 1). Then, taking into account relation ( |A.5[) , one can represent the 
Hamiltonian as 

H a = N + 1 + uU + . 

Therefore, the energy levels corresponding to the odd eigenstates \2l + 1) now are stable, 
whereas the shift of the levels corresponding to the even eigenstates \2l) depends on v. Special 
cases correspond, again, to the integer values of the deformation parameter, when RDHA 
specifies parabosons. For the parabosons of odd order p — v + 1 = 2k + 1, k — 1, 2, . . ., the 
spectrum is like in the non-deformed case v — but with the only difference: it contains k 
'holes' in its lower part in comparison with the bosonic system. For odd values v = 2k + 1, 
k = 1, . . ., the spectrum looks like the spectrum of usual supersymmetry (k = 1) shifted for 
AEi = E = 2, or like that of the nonstandard supersymmetry (k = 2,3, . . .). 

Let us consider the case of parabosons of even order p — v + 1 = 4, 6, . . ., and shift the 
Hamiltonian for the constant corresponding to the lowest energy level, H a — ► H® 



1 a ' 



Hi = a~a + - 2. (3.3) 

Then, in the case p = 4 the spectrum of the system looks like the spectrum of the superoscil- 
lator [4~2"| being the system in the phase of unbroken supersymmetry with the lowest singlet 
zero energy level and all other levels to be supersymmetric doublets. But what are the cor- 
responding conserved (commuting with H®) odd supergenerators and what superalgebra do 
they form? Such mutually conjugate generators for the general case u = p — l = 2fc + l are 
given by 

Q+ = (a + ) 2fc+1 n + , Q- = (a-) 2/£+1 IT_. (3.4) 
They satisfy the following relations: 

Q + Q- = n k (H° a )u_, q-q + = n k (H° a )U + , 

where 

H k (H° a ) = (H° a - 2k + 2) ■ (H° a - 2k + 4) . . . H° a . . . (H° a + 2k) ■ (H° a + 2k + 2). 

Using them, we arrive at the following nonlinear superalgebra which characterizes the system 
( p.3| ) at v = 2k + 1, k = 1,2,..., (or the system given by the Hamiltonian ( |3.1| ) with the 
substitution H° = H a - 2): 

Q ±2 = 0, {Q + ,Q-} = H k (H° a ), [fl2,g ± ] = 0. (3.5) 



On the other hand, with the 'first order' odd generators given by Eq. ( |3.2| ) the system fl3.3| ) 
can be characterized by the following central extension of the N — 1 superalgebra with the 
generators Q ± , H® and R: 

Q ±2 = 0, {R,Q ± } = 0, [R,H° a ]=0, 
{Q + ,Q-} = H° a -kR + k + 2, [H° a ,Q ± } = ±2kQ ± , [R,Q ± ]=±2Q ± , R 2 = 1. 

Therefore, in the case p = 4 the system (|3.3|) has the usual spectrum of the system in 
the phase of exact supersymmetry, but its odd conserving supergenerators form with the 
Hamiltonian the nonlinear superalgebra: 

{Q + ,Q-} = H° a (H a +2)(H° a +A), = 0, Q ±2 = 0. (3.6) 
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Note that here the ground (supersymmetric vacuum) state is \vac) = |1), and, so, it is odd, 
R\l) = — 11), whereas |0) belongs to the first excited supersymmetric doublet. Therefore, 
this system is not included in the minimally bosonized version of supersymmetric quantum 
mechanics where the supersymmetric ground state is always even ^(| and this explains 
the appearance of the nonlinear superalgebra. The systems with the Hamiltonian ( |3.3| ) 
corresponding to p — 6, 8, . . . are characterized by the same property: their ground state is 
odd: \vac) = |1), H®\vac) = 0. 

The coordinate realization ( |1.4|) , (|1.5|) , results in the following form for the Hamiltonian 
H® = a~a + — 2 and supercharges (|3.4|): 

h ° - H-^ + * 2+ ^- 3+ "(i +i H' (37) 

with v — 2k + 1. Therefore, among the Calogero-like 2-particle systems with exchange in- 
teraction given by the Hamiltonian ( |3.7D the cases corresponding to the odd values of the 
parameter v = 3,5,7,... are characterized by the presence of the hidden nonlinear super- 
symmetry, whereas the system given by the Hamiltonian H® + 2 at v — 1 possesses linear 
supersymmetry ( |1.3| ) in the spontaneously broken phase. In particular, as we observed, the 
Calogero-like system (|3.7|) corresponding to v = 3 is a rather peculiar one: its spectrum coin- 
cides exactly with the equidistant spectrum of usual linear superoscillator, but the associated 
supercharges ( |3. 8| ) satisfy in this case the nonlinear superalgebra fl3.6|). 



4 Realization of parabosons by two identical fermions 

Suppose that we have the system of two identical fermions such that its space motion is 
confined to the line whereas their spin degrees of freedom are not restricted and are given 
by the total vector spin operator with components Jj = |(<7j ® 1 + 1 ® cr^). Then, taking into 
account the Pauli principle and omitting the dependence on the center of mass coordinate, 
X = \{xi + X?), we describe the two-fermion system by the wave functions of the form [Mj 

*{x) = xi 3 *P J -(x) + XaMx). (4.1) 

Here j'3 = —1,0, 1, x — X\ — x-i is the relative coordinate; xt 1 = X^ 1 = l — )l — ) an d 

X° = "^(|+)| — ) + | — )|+)) are symmetric spin states forming a vector triplet, JiJiXl 3 — 2Xs 3 > 
JzX^s — J3Xs 3 5 whereas Xa = ^(1+) l — ) — l — i s antisymmetric spin-0 singlet state, 
JiXa = 0; ip- are °dd functions, ip^(-x) = — ip^(x), and ip + is an even function, ip+(— x) = 



ifj+(x). The physical operators are those transforming the states of the form fl4.1|) into the 
states of the same form. If we assume, moreover, that there is a rotational Jf -symmetry in 
the spin space, then the algebra of physical operators may be generated by the operators 

H 

A + = U U Pi + , A- = /_ (x, pi o_, 
\ ax \ ax 1 
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where f ± {-x, -£) = ±f±{x±); + = l,E 3 ,Ei, i = 1,2,3; G_ = E 1; E 



2- 



S l = ^(°3 ® 1 - 1 ® (^3) j S 2 = -((Ti ® 2 - CT 2 ® 0"l), S3 = ^{Vl ® &1 + 0"2 (g> 2 ) 



01 - 02 ® a 2 J 



—2 



-(0i 
2 V 1 



02 + 02 



Operators Ej and 

S 3 Xs(a) = +(-)Xs(a), S i X =t 
As 1 As 



501; 

X&; 



"3 



03 + 03 ® 1; 



/2 77. 



satisfy the relations EjEj = XJy + ie^E^, X = 1 
= 0. They act on the spin states in the following way: EiXs(a) — Xa(s), 



0, H lX ± 


As * 



X T , H 3 x = ±X , ^iXs(a) = 0, where we have used 

the notation x d 

Let us return to the i?-deformed Heisenberg algebra. In the coordinate representation an 
arbitrary function may be decomposed into even and odd parts, ^/(x) = ^/ + (x) + ^/_(x), 

map the even and odd functions into the 



x 



R*ff±(x) = ±ty±(x). Even operators / 4 

functions of the same parity nature, whereas odd operators f-(x, -|-) change the parity of 
functions. Then, taking into account the properties of operators Ej and E, specified above, 
one can arrive at the following identification of the states and operators of RDHA on the one 
hand, and, on the other hand, of the states and nontrivial physical operators in the js = 
sector of the system of two identical fermions: 





Fermion system (j 3 = 0) 


RDHA 


States 


Xs^- 






Xa1p+ 




Operators 




/+ 




;+s 3 






/-Ei 


f- 




/-E 2 


if-R 



This identification means that the operators 

1 



.4 



± 



V2 



Si 



x =F 




(4.2) 



together with the operator E 3 are the analogs of the creation-annihilation operators ( |L5l 
and of the reflection operator R. Indeed, the operators A ± and E 3 satisfy the relations 



[A-,A< 



;i + z/e 3 )x, {z 3 ,a ± } = o, 



y2 
^3 



X, 



which in the sector j'3 = are reduced to the relations of RDHA (|1.6| ). Therefore, all 
the obtained results on the supersymmetries of single-mode parabosonic systems can be 
translated to the js = sector of the system of two identical fermions. Then the Hamiltonian 



Hi = A + A- 



d 2 



dx 2 



W 2 +v-l-(W' + i/-l)E 8 



X 



(4.3) 



with W +>v 

described by the Hamiltonian 



+x — corresponds to Eq. (|2.2j) . So, the system of two identical fermions 

reveals linear and nonlinear supersymmetries in J3 = 
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sector when the parameter v takes the values v = 1 and v = 3, 5, 7, . . ., respectively. Note 
that here the states with j 3 = ±1 always have zero energy. Can we modify the Hamiltonian 
( |4.3| ) in such a way that the infinitely degenerated zero energy levels corresponding to js = ±1 
states would be unfrozen without changing the supersymmetric properties of j$ = sector? 
A natural modification of the Hamiltonian (|4.3|) is 



n 2 



d 2 



dx 2 



a* 



(4.4) 



The difference A = H^ — H£ is given by the operator 



A = - 

2 



dx 2 



+ W 2 +u-l-(W'+u-l) 



■ J 2 



which acts nontrivially only in the j% = ±1 sector. Moreover, this difference is chosen in 
such a way that its action on the states xf 1 ^- 1 ^)^ XJ'V- 1 ^) coincides with the action of 
the operator H£ on the state xfy-fe)- As a result, the energy levels of Hamiltonian ( |4.4|) 
corresponding to all odd eigenstates (given by odd functions of x) have a multiplicity equal 
to 3. The spectrum of the corresponding two-fermion system is shown in Fig. 2a for the case 
of v = 3. 

The two-fermion Hamiltonian 



Hi = A-A^ 



d 2 



dx' 



+ wl >v -v + i-(wL,„-i' + i)x s 



1 



(4.5) 



with = —x — j- is the analog of the Hamiltonian ( |3.1| ), and at u = 2k + 1 it reveals 

supersymmetry in j 3 = sector. The modified (and shifted for —2) Hamiltonian 



H°J 



1 



d 2 



dx 2 



+ W1 



1 



{WL v -v + l)-<Ti®ai (4.6) 



is different from Hi in 



A 



dx' 



+ W 2 v -u + l-{W'_ v -u + l) 



■Ji-2. 



The energy levels corresponding to odd eigenstates of the Hamiltonian ( |4.6| ) (represented by 
odd functions of x) have also the multiplicity equal to 3 and its spectrum is shown in Fig. 
2b for the case v = 5. 

Thus, the supersymmetries of the single-mode parabosonic systems may be realized in 
the sector j 3 = of the system of two identical fermions. The modified Hamiltonians ( [4.4| ) 
and (|4.6|) have the form of the Hamiltonians corresponding to the simplest 2-particle case 
of spin chain models with inverse-square interaction |^TJ. These Hamiltonians unfreeze the 
energy levels of the states with j 3 = ±1 (which for the Hamiltonians ( |4.3j ) and ( |4.5| ) turn out 
to be infinitely degenerated states of zero energy) leading to the multiplicity equal to 4 for 
higher degenerated states. 
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33 

a) 



±1 



3s 



33 

b) 



±1 



Figure 2: a) Spectrum of the two-fermion system ( |4.4| ) at v — 3; b) Spectrum of the two- 
fermion system (|4.6|) at v = 5. 

5 Nonlinear supersymmetry in boson-fermion system 

Nonlinear supersymmetry specific to the pure parabosonic system may also be realized by 
the boson-fermion system. To construct such a boson-fermion system with nonlinear super- 
symmetry, it is convenient to start from the (pseudo) classical level using the ideas of ref. 



45]. So, let us consider a nonrelativistic particle of unit mass carrying the odd degrees of 



freedom described by two Grassmann variables 6 a , a = 1, 2. The general Lagrangian of such 
a system is 

L = l -x 2 - V(x) - L(x)T + l -ej a , 
where V(x) and L(x) are two arbitrary functions and T = —i9i9 2 = 6 + 6~ , 9^ 



72^ 



h±i6 2 



The nontrivial Poisson-Dirac brackets for the system are {x,p} = 1, {0 a ,8t,} = —i8 a bi an -d 
the Hamiltonian H = |p 2 + V(x) + L(x)T generates the following equations of motion: 

x = p, p = -V'(x)-L'(x)r, 6 ± = ±iL(x)9 ± . 
Let us construct the oscillator-like variables 



with some function W(x). We find that 



— (W(x)Tip) 



±iW'B ± ±—(V'- WW) ± -=L'Y. 



V2 



V2 



As a consequence, the odd quantities 



<3fc = (B±) k 9T, 



where k is some natural number, satisfy the following equations of motion: 

Qt = ±iQt(kW -L)± -^=k(V - WW')ff*. 



(5.1) 
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They will be odd integrals of motion if we choose V = \{W 2 )' , and put L = kW. Therefore, 
when the functions V(x) and L(x) are related as L(x) = kW(x), V(x) = \W 2 {x) + C , where 
C is a real constant, then the odd integrals ( |5.1| ) are the integrals of motion additional to H 
and T. Together with the Hamiltonian 

Hk — -jip 2 + w 2 ) + kW'F + C (5.2) 

they form the following nonlinear superalgebra: 

{Qt Qk) = ~i(Hk - C)\ {H k , Q±} = {Q+, Qt} = {Q k ,Q k } = 0, (5.3) 

whereas T is the classical analog of the grading operator, {T,Q k } = {T,H k } = 0. So, the 
class of (pseudo) classical systems revealing nonlinear super symmetry algebra is given by an 
arbitrary function W{x): 

L = 1 (x 2 - W 2 (x) + ikW'(x)e ab 6 a 6 b + t6 a 6 a ) , (5.4) 

where we have put C = 0. 

Let us consider the simplest example corresponding to the case of linear superpotential 
W = x. The quantum Hamiltonian of such a system is the sum of the Hamiltonians of 
harmonic bosonic and fermionic oscillators with the fermionic frequency being the integer 
multiple of the bosonic one: 

H k = H b + kH f , H b = b + b-, H f = f + f-. (5.5) 

In the case k — 1 the system is a superoscillator for which the quantum analog of 
superalgebra ( |5.3| ) is given by Eq. (|1.3|). For k = 2,3, . . . the quantum analog of superalgebra 
( |5.3|) for the system ( |5.5|) takes the form similar to the nonlinear superalgebra of parabosonic 
system Q or (£3J): 

{Q+,Q k -} = H(H-l)...(H-k + l). 

The boson-fermion system (|5.5| ) possessing the nonlinear supersymmetry of the type we have 
revealed in the single-mode parafermionic systems is analogous to the system of bosonic 
oscillators with rational ratio of frequencies, which, in turn, is characterized by the nonlinear 
finite iy-symmetry 0, 0. 

For the superpotential W(x) different from the linear one the commutator of operators 
B + and B~ is nontrivial and in general case this leads to the quantum anomalies destroying 
the classical nonlinear supersymmetry |f^B| . 



The 'hole' structure specific for the parabosonic systems (|2.2|) and at v = 2k may 
also be reproduced by the simplest boson-fermion system in the obvious way: the spectrum 
of the Hamiltonian H = 2[b + b~ + {k + \)f + f~] reveals k 'holes' in comparison with the 
spectrum of the bosonic number operator H b = b + b~ . 

6 Discussion and outlook 

To conclude, let us discuss the obtained results and indicate the problems that deserve further 
attention. 
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We have found that the single-mode systems of parabosons of even order with the Hamil- 
tonians H n = a + a~ and H a = a~a + possess the usual or nonstandard supersymmetries. 
These supersymmetries are characterized by the superalgebra of the form ( |1.3| ) linear in the 



Hamiltonian, or by the nonlinear superalgebras (|2.7|) or (|3.5| ). The peculiar nature of the 
revealed supersymmetries is encoded in the structure of the corresponding conserved odd op- 
erators ( |2.6|) and ( |3.4|) . They contain the projectors n_i- , which in terms of the parabosonic 
creation-annihilation operators may be represented as 

IT + = cos 2 F, II_ = sin 2 F, F=j{a + ,a~}, 



where we have used Eqs. ( A. 3 ) and ( A.4 ) and taken into account that v — 2fc + l. Therefore, 
the conserved odd supercharges are essentially nonlinear operators of a ± . 

The i?-deformed Heisenberg algebra, with the help of which we have analyzed hidden 
supersymmetries of the single-mode parabosonic systems, can be characterized in accordance 
with Eqs. (|A.5| ), ( |A.4j ) by the relations 



a + a- = F(N), a~a + = F(N + 1), (6.1) 

where F(N) = N(—i) N + z/sin 2 The discussed supersymmetric structure of the para- 
bosonic Hamiltonians H n and H a is hidden in the relation 

F(2n + 1) = F(2n + v + l), v = 2k + 1, k — 0,1, (6.2) 

RDHA can be treated as a particular case of the so called generalized deformed oscillator 
algebra |43| given by the relations (|6.1| ), by commutation relations characterizing the number 
operator, [N, a^] = ±a ± , and by the structure function -F(O) = 0, F{n) > 0, n — 1, being 
an analytic function. Therefore, it is interesting to investigate in the context of hidden 
supersymmetries the general class of single-mode generalized deformed oscillator systems 
whose characteristic functions have the property of degeneracy analogous to ( |6.2j) . The first 
results in this direction have been obtained recently by observation of the hidden linear and 
nonlinear supersymmetries in parafermionic type systems (whose characteristic functions are 
characterized by additional relation F{p + 1) = 0) |47[ . 



We saw that the realization of the creation-annihilation parabosonic operators via Eqs. 
L4|) , ( |1.5|) results in representation of the Hamiltonians H n = a + a~ and H® = a~a + — 2 



in the form ( |2.8| ), fl3.7| ), which may be treated as 2-particle Calogero-like models with ex- 
change interaction. The obtained results say that the odd values of the "coupling constant" 
v characterize the special subclass of such Calogero-like models possessing the hidden su- 
persymmetries. These supersymmetries are described generally by nonlinear superalgebras. 
To our knowledge, the nonlinear supersymmetries is a new feature of Calogero-like systems 
with exchange interaction not observed earlier. Therefore, it is interesting to generalize the 
systems (|2.8|) , ( |3.7|) for the n-particle case, n > 2, (by introducing the set of corresponding 



exchange operators [|30], f|]) to reveal the possible hidden nonlinear supersymmetric structure 
in n-particle Calogero-like models with exchange interaction. As in the 2-particle case (|2.8p , 
(|3.7D, the peculiar nature of such supersymmetries has to be encoded in the 'bosonized' struc- 
ture of the supercharge operators which have to depend explicitly on exchange operators and 
generally will be higher order operators in derivatives d/dxi (see Eqs. (|2.9|), (|3.8D). 
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We have observed that the spectrum of supersymmetric parabosonic systems character- 
ized by the presence of several supersymmetric singlet states may be reproduced by the 
reduced parasupersymmetry. The interesting problem is finding the classical analog for the 
reduced parasupersymmetry. The corresponding classical systems may be realized with the 
usage of paraGrassmann variables, and the reduction conditions ( |2.4| ) will appear as the 
quantum analogs of classical first-class constraints. Such a scheme of obtaining nonlinear 
supersymmetries by reduction of the appropriate parasupersymmetric systems will be anal- 
ogous to the Poisson reduction of the Kirillov Poisson structure, with the help of which the 



nonlinear finite M^-algebras are constructed (TPI . 

We have shown that the .R-deformed Heisenberg algebra and the associated supersym- 
metries of the single-mode parabosonic systems are realizable in the j'3 = sector of the 
system of two identical fermions. The correspondence is one-to-one and the crucial role is 
played here by the Pauli principle. On the other hand, by analogy with relations ( f4.2|) , one 
can realize the operators and R as 



+ 



V2 



x =F 



v 

2x 



0-3 



R = a 3 . 



(6.3) 



Operators (|6.3|) satisfy the relations of RDHA ( |1.6| ), but such a realization is reducible. 
This is clear just from the observation that the operator 1Z = (— 1)^03, with N constructed 
according to Eq. ( |A.3| ), commutes with all the operators (|6~3|). 

Though we have seen that the nonlinear superalgebra may be reproduced via the appro- 
priate modification of the classical analog of the usual supersymmetric quantum mechanics, 
nevertheless the quantum systems obtained after quantization are different from the su- 
persymmetric pure parabosonic systems. Indeed, in supersymmetric quantum mechanical 
systems there are independent fermionic f ± and bosonic b ± operators. Proceeding from 
the parabosonic operators a^ 1 one can also construct the operators b and / satisfying the 
(anti) commutation relations [b~,b + ] = 1, {f~,f + } = 1, f ±2 = 0, but they will not be inde- 
pendent: [b^J^] ^ (see refs. |39|, |34|). One can show j|6| that the parabosonic systems 
with nonlinear supersymmetry may be related to the nonlinear supersymmetric quantum 



mechanical analog of classical systems ( |5.2| ) in the same way as the minimally bosonized 
version of supersymmetric quantum mechanics with its linear supersymmetry ( |1.3| ) is re- 
lated to the usual supersymmetric quantum mechanics, i.e. via the special nonlocal unitary 
transformation applied to the latter with subsequent reduction to a subspace of one of the 



eigenstates of the total reflection operator [40 



Among the quantum systems corresponding to the modified classical analog of super- 
symmetric quantum mechanics fl5.4Q , one can find anomaly- free subclasses of the systems 
with nonlinear superpotential [f|(J. Investigating them, we hope to establish a relationship 
between nonlinear supersymmetry observed here and the known construction of the hierar- 
chy of Hamiltonians related by linear supersymmetry, which plays important role in some 
integrable systems [[J . Since the nonlinear supersymmetry of the boson-fermion system ( |5.5p 
is a super-analog of the nonlinear symmetry of the system of bosonic oscillators with rational 
ratio of frequencies |2D|, [B5J, investigation of the anomaly-free quantum systems correspond- 
ing to the classical systems ( |5.4| ) with nonlinear superpotential W (2) attracts also attention 
from the view-point of nonlinear finite iy-symmetries. 
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A Appendix: Irreducible representations of RDHA 

Here we give some information on RDHA and its representations, which is necessary for 



understanding the main text. The details and further references may be found in E5 



As a consequence of Eqs. (E6), the operators a obey the trilinear commutation relations 



[{a~,a + },a ± \ = ±2a d 



(A.1) 



If one defines the vacuum state as a |0) = 0, R\0) = |0), (0 1 0) = 1, than Eq. (|L6|) gives the 
relation 

a - a + \0) = (l + u)\0). (A.2) 

For integer values of the deformation parameter v — p — 1, p — 1, 2, . . ., the relation ( |A.2| ) 
together with the trilinear relations ( |A1|) specifies the creation-annihilation operators of 
the parabosons of order p, whereas for non-integer v > — 1 (see below) the operators a 
satisfying relations (|1.6|) may be treated as the creation-annihilation operators of generalized 



parabosons [23 



The trilinear relations ( |A.1| ) can be represented equivalently as {G, a } = 0, G = 
[of, a + ] — 1. As a consequence, [G 2 , a ] = 0, and one concludes that in irreducible represen- 
tation G 2 = const. If, moreover, G is hermitian operator, then G = vR with {R,a ± } = 0, 
R 2 = 1 and v G R. In other words, the trilinear relations themselves lead to RDHA ( |1.6|) . 

Relations ( |A.1| ) mean that up to the additive constant the number operator is given by 
a half of anticommutator of a + and a~ . The commutator from Eq. (|1.6|) together with the 
specified above definition of vacuum state fixes the additive constant and we get the number 
operator 

N = \{a\a-}- l -{l + v). (A.3) 
In its terms, the reflection operator can be realized as 



R 



COS irN. 



(A.4) 



With relations (1^) and (|A.3|) , the normal and antinormal products of the creation- 
annihilation operators can be given in terms of the number operator as follows: 



a a 



a a 



N + 1 + zTL 



where U± = |(1 ± R) are the projector operators, U 2 . = n_, Hj_ = n + , n + n_ 



(A.5) 
= 0, 



L± - 2 , 

n+ + n_ = i. 

From the requirement (n\n) > 0, |n) oc (a + ) n |0), n = 0, 1, 2, . . ., it follows that RDHA has 
unitary Fock-type representations when u > — 1. In this case the orthonormal eigenstates 



of the number operator are given by \n) = C n (a + ) n \0), C n = ([n\ u \) 1 ' 2 , [ 



n 
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[l] u = I + |(1 — (—l) l )v. For integer values v = k, k = 0, 1, . . ., these representations are 
equivalent to the Fock representation for parabosons of the corresponding order p = k + 1. 

Let us turn to the Schrodinger (coordinate) representation. Here, in correspondence with 
Eq. (|1.5|) , we have \/2a~ = x + 4 — j-R. Then the ground state, specified by the conditions 
R^ = ^o, a~^o = 0, is given by the function ^o(x) = \x\ u ^ 2 exp(— x 2 /2)J\f, where J\f is a 
normalization constant. With the scalar product defined by 

/+oo 
Vl(x)V 2 (x)dx, (A.6) 



we get 



POO 

(tf # ) = 27V 2 / Id" exp(-x 2 )dx. 
Jo 



The integral is converging (at lower limit) iff v > —1, i.e. the ground state is given by 
the normalizable function only for v > — 1. Therefore, one concludes that in Fock and 
Schrodinger representations singularity of v < — 1 reveals itself in different ways. In Fock 
representation for v — — 1 the state |1) has a zero norm, (1|1) = 0, whereas for v < 
— 1 there appear the states with negative norm (in particular, (1|1) < 0). In Schrodinger 
representation, the ground state is not normalizable for v < — 1 with respect to the usual 
scalar product ( |A.6|) . 



For the sake of completeness we note that the algebra ( |1.6| ) reveals some universality 



25 , |35| : for the special values v = — (2p + 1), p — 1,2,..., it has non-unitary finite- 
dimensional representations characterized by the relations (a ± ) 2p+l = 0. One can pass over 
to the hermitian conjugate operators f~=a~,f + = Ra + . They generate the parafermionic 
algebra of order p = 2 (for v = —3), or supply us with the deformation of parafermions of 
order p — 4, 6, . . ., which contains the internal ^-grading structure: 

[[f + J-]f ± } = 2(2I 3T l)Rf ± , 

where 

r = (-1)^, j 3 = /-] . ( _i)§[/ + J-]+ P . 
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